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Abstract 



Tracer diffusion in a granular gas in simple shear flow is analyzed. The anal- 
ysis is made from a perturbation solution of the Boltzmann kinetic equation 
through first order in the gradient of the mole fraction of tracer particles. The 
reference state (zeroth-order approximation) corresponds to a Sonine solution 
of the Boltzmann equation, which holds for arbitrary values of the restitu- 
tion coefficients. Due to the anisotropy induced in the system by the shear 
flow, the mass flux defines a diffusion tensor Dij instead of a scalar diffusion 
coefficient. The elements of this tensor are given in terms of the restitution 
coefficients and mass and size ratios. The dependence of the diffusion tensor 
on the parameters of the problem is illustrated in the three-dimensional case. 
The results show that the influence of dissipation on the elements is in 
general quite important, even for moderate values of the restitution coeffi- 
cients. In the case of self-diffusion (mechanically equivalent particles), the 
trends observed in recent molecular dynamics simulations are similar to those 
obtained here from the Boltzmann kinetic theory. 

PACS number(s): 05.20.Dd, 45.70.Mg, 51.10.+y, 47.50.+d 
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I. INTRODUCTION 



Granular systems under rapid flow conditions can be modeled as a fluid of inelastic hard 
spheres. In the simplest model the grains are taken to be smooth so that the inelasticity is 
characterized through a constant coefficient of normal restitution. The essential difference 
with respect to molecular fluids is the absence of energy conservation yielding modifications 
of the usual hydrodynamic equations. Due to the kinetic energy dissipation in collisions, 
energy must be externally injected to the granular gas in order to achieve a stationary state. 
In some experimental situations, the granular system is driven into flow by the presence of 
a shear field. In this case, a steady state is possible when the amount of energy supplied 
by shearing work is balanced by that lost due to the inelastic cooling. The study of the 
rheological properties of this steady shear flow state has received a great attention in the 
past years, especially in the case of one- component systems However, much less is known 
in the more complicated case of multicomponent mixtures of grains. 

An interesting problem is the analysis of diffusion in granular shear flows. The under- 
standing of mass transport in granular systems is of practical interest since, for instance, 
powders must frequently be mixed togheter before any sort of processing can begin. The 
self- diffusion phenomenon in granular flows has been studied earlier. Experimental studies 
include both systems with macroscopic flows [^,0 and vertically vibrated systems 0. Com- 
plementary computer simulation studies have been also carried out ||, with special emphasis 
on the influence of the solid volume fraction on the diffusive motion of the grains. In general, 
all these previous studies were limited to observing the diffusion in only one direction, usu- 
ally the direction parallel to the velocity gradient. However, due to the anisotropy induced 
in the system by the presence of shear flow, a diffusion tensor is required to describe the dif- 
fusion process instead of a single diffusion coefficient. To the best of my knowledge, the only 
attempt to measure the elements of this self-diffusion tensor has been made by Campbell 
||. He measured these elements via molecular dynamics simulations by using both particle 
tracking and through velocity correlations. Both methods were found to agree with rea- 
sonable accuracy. Very recently JF|, the self-diffusion coefficients have been experimentally 
measured in a granular system under Couette flow by employing image technology. 

In the context of kinetic theory, studies on granular flows in mixtures are scarcer. Most 
of them H are based on a Navier-Stokes description and, therefore, they are restricted to 
small velocity gradients, which for the steady simple shear flow is equivalent to the low- 
dissipation limit. For this reason, the diffusion is only characterized by a single coefficient 
which is not affected by the presence of the shear field. In addition, although these studies 
permit in principle different temperatures for the two species, they assume equal partial 
granular temperatures Tj in the quasielastic limit. Nevertheless, given the intrinsic con- 
nection between the shear rate and dissipation in this problem, energy nonequipartition is 
expected as the restitution coefficient decreases. As a matter of fact, some recent results 
obtained in molecular dynamics simulations of granular sheared mixtures || as well as in 
real experiments of vibrated mixtures in three jlTJ and two dimensions [TL\ clearly show the 
breakdown of energy equipartition. This implies that the temperatures T are different from 
the mixture temperature T. The consequences of this effect on the transport properties are 
in general significant, as has been recently found in the freely cooling case |]T~2| , p~3|| . In conclu- 
sion, a consistent theory describing diffusion in granular shear flows must take into account 
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both the tensorial character of the mass transport as well as the possibility of temperature 
differences. 

The aim of this paper is to get the diffusion tensor in a binary granular mixture under 
simple shear flow in the framework of the Boltzmann equation. Due to the complexity of 
the general problem, here we consider the special case in which one of the components (say 
for instance, the species 1) is present in tracer concentration. The tracer problem is more 
amenable to analytically treatment. First, the tracer particles are directly enslaved to the 
granular gas and, second there are fewer parameters [|T4| p76|] . Therefore, in this situation one 
can assume that the velocity distribution function f 2 of the excess component (granular gas) 
obeys a (closed) nonlinear Boltzmann equation while the velocity distribution function f\ of 
the tracer particles (impurities) satisfies a Boltzmann-Lorentz equation. The starting point 
of our study is a recent solution of the set of Boltzmann equations for a binary mixture of 
inelastic hard spheres under shear flow | 17fl . The corresponding Boltzmann-Lorentz equation 
for the impurities is solved by means of a perturbative scheme in powers of the gradient of 
the mole fraction of tracer particles. The main feature of this expansion is that the reference 
state around which we perturb is not restricted to small values of the shear rate, which for 
the steady shear flow problem is equivalent to arbitrary degree of dissipation. In the first 
order of the expansion, the tracer diffusion tensor is identified from the mass flux. Explicit 
expressions for the nonzero elements of this tensor are obtained by using a first Sonine 
polynomial approximation. These elements are given in terms of the restitution coefficients 
and the parametes of the mixture (masses and sizes) . 

The plan of the paper is as follows. In Sec. || we introduce the set of coupled Boltzmann 
equations describing the mixture and state the problem we are interested in. The state of 
the mixture in the absence of diffusion is analyzed and in particular, the nonzero elements of 
the pressure tensor P 2 of the granular gas are obtained in the leading Sonine approximation. 
Comparison of these results with previous theories is also presented. The Section ends 
studying the state of tracer particles with special emphasis in the evaluation of the specific 
dependence of the temperature ratio T\/T 2 on restitution coefficients, mass ratio, and size 
ratio. Section |TT| deals with the perturbation scheme used to solve the Boltzmann-Lorentz 
equation of the tracer particles when the diffusion takes place in the system. This Section 
contains the main results of the paper since we determine the tracer diffusion tensor in the 
first order of the expansion of the concentration gradient. The dependence of the nonzero 
elements of this tensor on the different parameters of the problem is illustrated in the three 
dimensional case, showing a good qualitative agreement with Campbell's simulations ||. 
Finally, in Sec. [IV] we close the paper with some concluding remarks. 



II. DESCRIPTION OF THE PROBLEM: GRANULAR MIXTURE IN SIMPLE 

SHEAR FLOW 

We consider a granular binary mixture composed by smooth inelastic disks [d = 2) 
or spheres (d = 3) of masses mi and m 2 and diameters <j\ and o 2 . Collisions between 
particles are inelastic and characterized by three constant restitution coefficients an, a 22 , 
and a\2 = a 2 \, where otij < 1 refers to the restitution coefficient for collisions between 
particles of species i and j. In the low-density regime, the one-particle velocity distribution 
functions /j(r, Vi;t) (i — 1,2) obey the set of nonlinear Boltzmann kinetic equations: 
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dt 



+ Vl • v) fi(r, v i; t) = £ Jij [vi|/iW, £(*)] 

/ 3 



(1) 



where the Boltzmann collision operator Jij [^i\fi, fj] describing the scattering of pairs of 
particles is 



x 



a<//i( r >vi,*)/.,-(r,v£,t) - /i(r,vi,t)/j(r,v 2 ,t) . (2) 



Here, d is the dimensionality of the system, <jy = (er, + dj) /2, is a unit vector along their 
line of centers, G is the Heaviside step function and gi 2 = Vi — v 2 . In addition, the primes 
on the velocities denote the initial values {vj,v 2 } that lead to {vi,v 2 } following a binary 
collision: 

vi = vi - fiji (l + o^ 1 ) (a- ■ gia)ff, v 2 = v 2 + pij (l + a^ 1 ) (ct • g 12 )a- , (3) 

where pij = rrii/ (rrii + rrij). At a hydrodynamic level, the relevant quantities are the number 
densities rij, the flow velocity u, and the "granular" temperature T. They are defined in 
terms of moments of the distribution ft as 



n; 



J dvfi(v), pu = J2pi u i = J dvrriivfi(v), (4) 



nT = Y,nm = Y,j rfv^\/ 2 /,(v), (5) 

i i 

where n = ri\ + n 2 is the total number density, p = p\ + p 2 is the total mass density, and 
V = v — u is the peculiar velocity Equations (H) and (|5|) also define the flow velocity Uj and 
the partial temperature Tj of species z, which measures the mean kinetic energy of species i. 

The collision operators conserve the particle number of each species and the total mo- 
mentum, but the total energy is not conserved. This implies that 

J2 J dY\m i V 2 J ij [Y\f h ft) = -fnTC , (6) 

where ( is identified as the "cooling rate" due to inelastic collisions among all species. At a 
kinetic level, it is also convenient to discuss energy transfer in terms of the "cooling rates" 
d for the partial temperatures Tj. They are defined as 

& = --^ E / dvmVj^Wi, ft) • (7) 
The total cooling rate ( can be expressed in terms of Q as 

C = r- 1 5> i :z&, (8) 

i 

where Xi = rti/n is the mole fraction of species i. 
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Let us now describe the problem we are interested in. We consider a granular binary 
mixture in which the masses and sizes of both species are arbitrary. Our aim is to analyze a 
diffusion problem in a mixture of grains driven into flow by the action of shear forces. Here, 
we consider the simple case of the tracer limit, namely, a binary mixture in which one of 
the components (say, for instance, 1) is present in tracer concentration (x\ < 1). In the 
tracer limit, one expects that the state of the granular gas is not affected by the presence 
of the tracer particles so that its velocity distribution function f 2 obeys a (closed) nonlinear 
Boltzmann equation. Furthermore, the mole fraction of tracer particles is so small that their 
mutual interactions can be neglected in the kinetic equation of J\. As a consequence, the 
velocity distribution function of tracer particles j\ satisfies a (linear) Boltzmann-Lorentz 
equation. Let us start describing the state of the mixture in the absence of diffusion. 



A. Granular gas (excess component) 

We assume that the granular gas is subjected to the simple shear flow. From a 
macroscopic point of view, this state is characterized by a constant linear velocity pro- 
file u = u 2 = a • r, where the elements of the tensor a are au = a5k x Se y , a being the constant 
shear rate. In addition, the partial density n ~ n 2 and the granular temperature T ~ T 2 are 
uniform. The temporal variation of the granular temperature T 2 can be obtained from the 
Boltzmann equation ([!]) as 

= ~ aP 2,xy ~ ^2?2, (9) 



where p 2 = n 2 T 2 , 



P 2 = m 2 J (ivVV/ 2 (10) 



is the pressure tensor of the gas, and in the tracer limit 

C 2 = --^jr J dvm 2 V 2 J 22 \v\f 2l f 2 ] . (11) 

The balance of energy (|9|) shows the different nature of this state for molecular and granular 
systems. While for elastic fluids (( 2 = 0) the temperature increases monotonically in time 
due to the viscous heating term aP 2}Xy , a steady state is possible for granular systems when 
the viscous heating is exactly compensated by the collisional cooling term (d/2)p 2 ( 2 . In that 
case, the shear stress P 2;Xy and the cooling rate ( 2 are related by 

aP2,xy = ~(2P2- (12) 

As a consequence, for a given shear rate a, the (steady) temperature T 2 is a function of the 
restitution coefficient a 22 . This steady state is what we want to analyze. 

The simple shear flow becomes spatially uniform when one refers the velocities of the 
particles to a frame moving with the flow velocity u: f 2 (r, v) — > / 2 (V). Consequently, the 
stationary Boltzmann equation for the excess component becomes 
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-aV y ^-f 2 (\) = J 22 [V\f 2 J 2 ]. (13) 

We are mainly interested in computing the nonzero elements of the pressure tensor P 2 . 
These elements can be obtained by multiplying the Boltzmann equation ( |T3"D by m 2 V k Vf 
and integrating over V. The result is 

m*- 2, km A ke , (14) 

where 

A M = m 2 J dVV k V e J 22 [V\f 2 , f 2 \. (15) 

To get an explicit expression for P 2 one needs to compute the right hand side of the set 
of equations (0). This requires the explicit knowledge of f 2 , which is not known even in 
the elastic limit. However, one expects to get a good estimate of the low moments of the 
Boltzmann collisional operator by expanding f 2 in Sonine polynomials and then truncate 
the series after the first few terms. This approach is similar to the usual moment method for 
solving the Boltzmann equation in the elastic case. Therefore, we take the leading Sonine 
approximation: 



/ 2 (V) - / 2 , M (v) i + H ( - s u j (v k v e - -^H kl 

where f 2 m is a Maxwellian distribution at the temperature of the gas, i.e 



(16) 



With the approximation (|16"D , the integrals appearing in the expressions of the cooling rate 
( 2 and the collisional moment A k e can be explicitly evaluated. The details of the calculation 
are given in the Appendix. This allows us to get the explicit expressions of the non-zero 
elements of P 2 . They are given by 

P * _ P * _ _ P * __ d+l + {d- l)a 22 

^, y y-^-----^,dd- 2d + 3-3a 22 ' [ 8J 

2 ' xy ~ (l + a 22 )(2d + 3-3a 22 y a ' [U) 

Pl xx = d-(d- l)Pl yy , (20) 

, 2 = d + 2 (1 + a 22 )(2d + 3 - 3a 22 ) 2 (l - a\ 2 ) 
a 32d d + 1 + (d - l)a 22 ' 1 ' 

Moreover, the (reduced) cooling rate Q * s 
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(22) 



Here, P2 = P2/P2, a* = a/v, Q = (2/v, and v is a characteristic collision frequency given by 
v = P2/1I1 with rj is the shear viscosity coefficient of the gas in the elastic limit, i.e., 
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i±l^-^Y{d/2)a 2 ^ l \m 2 T 2 yl\ (23) 
o 

Equations (|l8l )-(|2~2D generalize previous results derived in the three dimensional case 
|17|| . The expression fl2T|) clearly indicates the intrinsic connection between the velocity 
gradient and dissipation in the system. As a matter of fact, given that a 2 2 < 1, the range of 
(reduced) shear rates are defined in the interval < a* 2 < (d + 2) (3 + 2d) 2 /32d(d + 1). The 
parameter a* can be considered as the relevant nonequilibrium parameter of the system. In 
the elastic limit (0:22 — 1 which implies a* = 0), the equilibrium results of the molecular gas 
are recovered, i.e., P 2M = 8m- As said in the Introduction, the steady simple shear flow for a 
monocomponent granular fluid has been the subject of many previous works. Two interesting 
studies have been carried out by Jenkins and Richman for smooth inelastic disks |18| and 



by Brey et al. [19[ for a d dimensional system. The latter description has been subsequently 



extended to dense gases [^0] . The approximated theory of Jenkins and Richman |0| is based 
in a generalized Maxwellian distribution to model the steady shear flow. On the other hand, 
Brey et al. ]T9| solved a kinetic model equation of the Boltzmann equation and compared 
their predictions for the reduced temperature, shear stress, and normal stress differences 
with Monte Carlo simulations. Comparison between kinetic model results and simulation 



shows in general a good agreement. Our results differ from those reported in Refs. 0] and 
19| , although for practical purposes, the discrepancies between the different approaches are 
quite small, even for moderate values of «22- As an illustration, in Fig. |l] we compare the 
different theories for the (reduced) elements of the pressure tensor P2 in a two-dimensional 
system [d = 2). It is seen that the agreement is remarkable, although the discrepancies 
slightly increase as the restitution coefficient decreases. 



B. Tracer component 

In the absence of diffusion, the velocity distribution function fi of tracer particles satisfies 
the steady kinetic equation 

-^/i = W./i (24) 

The most relevant moment at this level of approximation is the pressure tensor Pi defined 

as 

P 1 =m 1 J rfVW/i(V), (25) 
The nonzero elements of Pi obey the set of equations: 

Q>kmP\,lm + 0>emPl,km = P>kl (26) 
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where 



B M = m 1 J dW k V t J 12 \y\h t f 2 }. 



(27) 



As done before in the case of P2, we estimate the collisional moment Bm by taking the 
leading Sonine approximation of f% 



/i(V) - /i,m(V) 



2Ti I pi 



1 

d 



v k v e - -y 2 s ke 



(28) 



where p\ = n{F\ and now jx t M is a Maxwellian distribution at the temperature of the tracer 
particles Ti, i.e., 



(29) 



As will be shown later, the partial temperatures T% and T2 are in general different so that 
the granular energy per particle is not equally distributed between both species. 

Once the collisional moment B^ is determined (see the Appendix), the (reduced) non- 
zero components of P* = P\/x\P2 can be easily obtained from Eq. (p6|). These components 
can be written in terms of the temperature ratio 7 = Ti/T 2 , the restitution coefficients a 22 
and a±2, and the parameters of the mixture. After some algebra, one gets 



p* p* 

r l,yy - 1, 



p* 



dd 



F + ^2, W 

G 



(30) 



p* 



a*^* - P 2 * 



G 



(31) 



where 



A% = d7 - (d - 



.v.y ' 



1/2 



# 3 



1 + M(d_l)(l + 0)(1 + Qtla ) 



(32) 



(33) 



G 



1/2 



^21 



;i + «i 2 ) 



_ V2 /^y- 1 
4^ U 2 J ' " v^c 1 + 

x {2[(d + 2)0 + d + 3] - 3/i 2 i(l + 0)(1 + a a2 )} 



(34) 



H = 1Z & 1 ^ (wTe)) (1 + ai2) [3/i2l(1 + m + ai2) - 2] • (35) 

Here, 9 = miT 2 /m 2 Ti is the mean square velocity of the gas particles relative to that of the 
tracer particles. To close the problem at this stage of approximation, it still remains to get 
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the temperature ratio 7. It can be obtained for instance, from the requirements 
corresponding counterpart for species 1. This yields 
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S.2- r l,. 



where the cooling rate Q 



Sl r 2,xy 

Qi/v for the tracer particles is (see the Appendix) 



(d + 2)V2 /oiA 



d-l 



1/2 



Ad 



V cr 2 / 



/i21 



(l + «i 2 ) [2-/i 21 (l + e)(l + a 12 )] 



and its 



(36) 



(37) 



The solution to Eq. ( |3B| ) gives 7 as a function of the restitution coefficients a 22 and a.12 
and the mechanical parameters of the mixture, i.e., the mass ratio fi = 1711/1712 and the 
size ratio w = 0\jo<i- Except for some limiting cases, Eq. (^) must be solved numerically. 
Thus, in the elastic case (a 22 = «i2 = 1), we recover the well-known equilibrium results with 
7 = 1 and 9 — > mi/1712 as required by the equipartition theorem. In the case of mechanically 
equivalent particles (mi = m 2 , a 22 = «i 2 , ci = o" 2 ), Eqs. (|1"8|)-(|22|) and (|30|) (p5|) lead 
to P2 = Pi and C 2 = C*> so that 7 = 1. Beyond the above limit expected, our 

results yield 7 7^ 1. The violation of energy equipartition in driven granular mixtures has 
been recently observed in molecular dynamics simulations of sheared mixtures || and in 
real experiments |TD|JTT|. This effect is generic of multicomponent granular systems and 
is consistent with previous results derived in the unforced case [p^ ,|l4]-|T6|l . To the best 
of my knowledge, the only previous theories including temperature differences have been 
proposed by Jenkins and Mancini Ell] and by Huilin et al. ||22||. However, both works are 



phenomenological with no attempt to solve the kinetic equation. Instead, they assume that 
the velocity distribution function is a local Maxwellian. This is reasonable for estimating 
the dense gas collisional transfer contributions to the fluxes, but not for evaluating their 
kinetic contributions. Both theories are applicable to a general flow field. In particular, the 
results obtained by Jenkins and Mancini for the temperature ratio in the low density 
limit for inelastic disks can be written as 
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1 - 



1 +n 
1 + w 



l + w 
2\x ' 



tti 2 J 



//i(l + /i) 



1 - a 22 ) 



(38) 



In Fig. we plot the temperature ratio 7 versus the restitution coefficient a for a size ratio 
w = 2 and three values of the mass ratio fi in the two-dimensional case (d = 2). For the 
sake of simplicity, henceforth we will assume that the spheres or disks are made of the same 
material, i.e., a = a 22 = «i 2 . Also for comparison, we show the prediction given by Eq. (|38"D 
in the case fi = 5. Important discrepancies between both theories appear even for values of 
a close to 1. As a matter of fact, the theory of Jenkins and Mancini predicts a violation of 
energy equipartition much more significant than our theory. It must be remarked that the 
quantitative predictions of our theory at the level of the temperature ratio have been recently 



confirmed by Monte Carlo simulations |T^,|23[. Regarding the influence of the parameters 
of the mixture, we observe that for large mass ratios the temperature differences are quite 
significant, even for moderate dissipation (say a ~ 0.9). The temperature of the tracer 
particles is larger than that of the excess species when the tracer grains are heavier than 
the grains of the gas. This behavior has been also found in the recent computer simulations 
carried out in rapid shear flow ||. 
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III. TRACER DIFFUSION UNDER SIMPLE SHEAR FLOW 



We want to study the diffusion of tracer particles immersed in a bath (granular gas) 
subjected to the simple shear flow. The diffusion process is induced in the system by a 
weak concentration gradient Vxi. However, given that the strength of the shear rate a is 
arbitrary, the mass flux (which is generated by the gradient Vxi) can be modified by the 
presence of the shear flow. As stated above, in the tracer limit the state of particles of species 
1 is mainly governed by the collisions with particles of species 2, so that the self-collisions 
among particles 1 can be neglected. Thus, the kinetic equation governing the evolution of 
the velocity distribution function f x reads 

di fl ~ aVy W h + {Vk + ahere) dr~ k fl = Jl2 [ v i /l > ^ ' (39) 

where here the derivative d/drk is taken at constant V. Tracer particles may freely exchange 
momentum and energy with the particles of the granular gas and, therefore, these are not 
invariants of the collision operator Jv}[fi-> f?\- Only the number density of tracer particles is 
conserved. More specifically, the mole fraction x% obeys the conservation law 

■j- + a ke r e —) Xl + = 0, (40) 
at or k J m 1 

where the mass flux ji is defined as 

ix = m x J dVVft(V). (41) 

When a 2 2 = «i2 — 1 (which gives a* = 0), the well-known Fick law establishes a linear 
relationship between the mass flux ji and the concentration gradient Vxi. This law defines 
the diffusion coefficient. For finite values of the (reduced) shear rate ajv (which means 
a 22 7^ 1), one expects that a generalized Ficks's law holds but now a diffusion tensor rather 
than a scalar should appear. Our aim is to get this tensor in terms of a 22 , a X2 , yU, and w . To 
this end, and assuming that the mole fraction X\ is slightly nonuniform, we solve Eq. ( |39|) 
by means of a perturbation expansion around a nonequilibrium state with arbitrary shear 
rate, which is equivalent to strong dissipation in the simple shear flow [see Eq. (]21f)1- Thus, 
we write 

fi = /i (0) + A (1) + ■ ■ ■ , (42) 

(k) 

where f x is of order k in Vxi but applies for arbitrary degree of dissipation since this 



distribution retains all the orders in a. The solution (f42|) qualifies as a normal solution since 
all the space and time dependence of fi occurs entirely through xi(r; t) and their gradients. 
The zeroth-order approximation corresponds to the simple shear flow distribution but 
taking into account now the local dependence on the mole fraction x\. Although the explicit 
form of fi°^ is not exactly known, only the knowledge of its second-degree moments (related 
to the pressure tensor P x ) is necessary to get the diffusion tensor in the Sonine approximation. 

The kinetic equation for f± can be obtained from the Boltzmann-Lorentz equation (|39D 
by collecting all the terms of first order in Vx x : 
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|A (0) - aV y ^-J? + ( Vi + a k ^jT = J 12 [V|A«, h ] . (43) 
According to the balance equation (£|(]), one nas 

d/i (0) = df[ 0) dx 1 = dndfP 
dt dx\ dt dr k dxi 



where use has been made of the fact that the zeroth-order approximation to the mass flux 

df® dtf) dx x 



vanishes, i.e., jj° = 0. Moreover 



drk dxi drk 
Using (pHj ) and (pE5|), Eq. (^3|) can be written as 



^+ a )-^¥ v - <49) 



(45) 



(^ + a)/1 1 »^(V.V, 1 ), (46) 

where A is the Boltzmann-Lorentz collision operator 

AA W = Ji 2 [V|A (1) ,/ 2 ]. (47) 
The solution to Eq. (|46| ) is proportional to Vxi, i.e., it has the form 

f[ 1 \v)=A(V)-Vx 1 (48) 
Substitution of this into Eq. ( [46"1) yields 



The first order approximation to the mass flux is given by 

j^=mi J dVVA (1) (V) 

= -D • Vxi, (50) 

where the tracer diffusion tensor is 

D M = -mi J dVV k A e (V). (51) 

The solution to the integral equation ( f4"9| ) allows one to determine the quantity A. From 
this solution one can determine the tracer diffusion tensor by means of Eq. (]51f) . 

In order to get an explicit expression for the tensor D we need to know the quantity A. A 
good estimate of A to evaluate the mass flux is given by the first Sonine approximation, 
in which only the leading term in the expansion of A(V) in Sonine polynomials is kept. 
Thus, we take the following approximation to A: 
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A(V) 



■V-D/ 1)M (V). 



Using Eq. (^), an equation for the diffusion tensor is easily derived from Eq. 
result can be written as 



(52) 
). The 

(53) 



where the nonzero components of P* are given by Eqs. (|30|) -(32) and we have introduced 
the tensorial quantity 



(54) 



The expresion of the tensor O has been obtained in the Appendix with the result 



° = \ ^l^^ia^T./m^il + a l2 ) [(1 + OW 1 ' 2 
The solution of Eq. ( |53"D is 

D = p 2 (a + Siy 1 ■ P* x 



d + 2 d + 2 2 



(55) 



(56) 



Equation (|56| ) is the primary result of this paper. It provides an explicit expression of the 
tracer diffusion tensor of a granular binary mixture in simple shear flow. The elements of 
this tensor give all the information on the physical mechanisms involved in the diffusion of 
tracer particles in a strongly sheared granular gas. In the absence of shear rate (what is 
equivalent to a 22 = o>\ 2 = 1), D ki = D 5ki, where 



D 



d r(d/2) 
4 v / 2vr( rf - 1 )/ 2 



d-i 

'12 



/i(l + /i) (m 2 T 2 



,1/2 



(57) 



is the tracer diffusion coefficient of a molecular gas [p4|j . As the restitution coefficients 
decrease, rheological effects become important and the elements of the diffusion tensor are 
different from the one obtained in the equilibrium case. The dependence of the diffusion 
coefficients on the restitution coefficients a 22 and a± 2 as well as on the mass ratio fi and 
the size ratio w is highly nonlinear. As happens for elastic fluids ||25| , Eq. (|56| ) shows that 
diffusion under simple shear flow is a very complex problem due basically to the anisotropy 
induced in the system by the shear flow. 

To illustrate the dependence of the elements D^e on the parameters of the problem, let 
us consider a mixture of inelastic hard-spheres (d = 3). According to Eq. (0), D xz = D zx = 
D yz = D zy = 0, in agreement with the symmetry of the problem. Consequently, there are five 
relevant elements: the three diagonal and two (D xy , D yx ) off-diagonal elements. In general, 
D xx 7^ D yy ytz f) zz and D xy ^ D yx . The off-diagonal elements measure cross effects in the 
diffusion of particles induced by the shear flow. Thus, for instance D xy gives the transport 
of mass along the direction of the flow of the system (x axis) due to a concentration gradient 
parallel to the gradient of the flow velocity (y axis). Both off-diagonal elements are negative. 
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Before analyzing the influence of the mechanical parameters of the mixture on diffusion, 
it is instructive to explore the particular case of self-diffusion, i.e., when the tracer parti- 
cles are mechanically equivalent to the gas particles. This situation involves only single- 
particle motion and it is therefore somewhat simpler to compute the diffusion coefficients. 
In particular, the temperature of the tracer particles is the same as that of the gas parti- 
cles and so 6 = 1 in Eqs. ©-© and (||). In Fig. |, we plot D* x - D* y , D* zz - D* y , 
(D xx + D* y + D* zz )/3 = (l/3)D* kk , -D xy , and — D* x as functions of the restitution coefficient 
a = «i 2 = a 22 . Here, DL = D^/Dq, with D given by Eq. fl57|) . We see that the deviation 
from the functional form for elastic collisions is quite important even for moderate dissipa- 
tion. Thus, for instance at a = 0.8, D xx - D* yy ~ 0.76, D* zz - D* yy ~ 0.046, (l/3)D* kk ~ 1.18, 
—D* ~ 1.039, and — D* ~ 0.42. The figure also shows that the anisotropy of the sys- 
tem, as measured by the differences D* x — D yy and D* zz — D* y , grows with the inelasticity. 
This anisotropy is much more important in the plane of shear flow (D* x — D* ) that in the 
plane perpendicular to the flow velocity (D* z — D* ). This is basically due to the fact that 

/ / 1 «-iti, . 1.2. 

As said before, Campbell |J has carried out molecular dynamics simulations to measure 
the nonzero elements of the self-diffusion tensor. In his work, the self- diffusion coefficients 
were nondimensionalized by the product of the shear rate and the particle diameter. In our 
units, this corresponds to the reduced tensor Djj = D*-/a*. Although the solid fractions 
analyzed in his simulations prevent us in general from making a quantitative comparison 
between our theory (restricted to dilute gases) and his computer simulations, we observe that 
the general qualitative dependence of the self-diffusion tensor on dissipation agrees well with 
our results, at least for the lowest solid fraction considered. Thus, theory and simulation 
predict that the magnitude of the normal diffusion coefficients follow the pattern D xx > 
D zz > D yy while, in general, the elements decrease as the restitution coefficient decrease. 
An exception to the latter rule is the element D xx , which does not depend sensitively on a. 
On the other hand, in Campbell's simulation work H, he found that the values of —D xy 
were roughly of the same magnitude as D yy provided that the solid fraction is smaller than 
0.4. This trend is not completely followed by our theory since the values of —D xy and 
Dy y are significantly different for highly inelastic spheres. Thus, for instance, at a = 0.8, 
-D xy ~ 2.43 and D yy ~ 2.15 but -D xy ~ 2.54 and D yy ~ 1.10 at a = 0.4. 

The dependence of the diffusion coefficients D*, on the restitution coefficient for different 
values of the mass ratio /i is illustrated in Fig. f|. In this case,we take a size ratio w = 2 
and two values of the mass ratio: /x = 2 and /i = 4. For a given value of the inelasticity, 
we observe that the deviations from the elastic results are more important as the tracer 
particles are heavier than the gas particles. 

IV. CONCLUDING REMARKS 

In this paper, we have described diffusion of tracer particles in a granular gas subjected 
to the simple shear flow. We have been interested in the steady state where the effect 
of viscosity is compensated for by the dissipation in collisions. Under these conditions, 
the resulting diffusion is anisotropic and, thus, cannot be described by a single diffusion 
coefficient. Instead, it must be described by a diffusion tensor whose explicit determination 
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has been the main objective of this work. In order to capture the essential aspects of such 
a nonlinear problem, we have considered here a granular mixture in the low density regime 
as a prototype granular system, which lends itself to a detailed description by means of the 
nonlinear Boltzmann kinetic equation. 

We have been concerned with the physical situation where a weak concentration gradi- 
ent coexists with a strong shear rate, which for the steady simple shear flow problem means 
strong dissipation. For this reason, the Boltzmann-Lorentz equation corresponding to the 
tracer particles has been solved by means of a perturbation expansion around a nonequi- 
librium sheared state. This implies that the different approximations of this expansion are 
nonlinear functions of the restitution coefficients as well as of the parameters of the mix- 
ture (mass and size ratios). To get explicit results, we have used a first Sonine polynomial 
approximation to evaluate the cooling rates and the collisional moments of the Boltzmann 
operators. The reliablity of this approximation has been recently assessed in the (pure) shear 
flow problem where it has been shown to agree very well with Monte Carlo simulations in 



the case of hard spheres [17 



The kinetic theory results show that the elements Dij of the diffusion tensor present 
a complex dependence on the restitution coefficients 022 and au and on the mass ratio 
\i = m\jm% and the size ratio w — 01/02. In the elastic case, = 1, = D 6ij, where 
D is given by Eq. (|57|) and one recovers the expression of the diffusion coefficient for normal 
fluids. The deviations of the tensor from the scalar D have two distinct origins. First, 
the presence of shear flow gives rise to the new tensorial term a + Q on the left hand side 
of Eq. (p3|) instead of the corresponding collision frequency of the elastic diffusion problem. 
Second, given that the tracer and fluid particles are mechanically different, the reference state 
(zeroth-order approximation of the expansion) of tracer particles is completely different from 
that of the gas particles. In particular, when /i 7^ 1 and/or w / 1, the temperature ratio 
7 is clearly different from 1 (as can be seen in Fig. ||]), confirming the breakdown of the 
energy equipartition. The effect of different temperatures for the tracer and gas particles is 
expressed by the appearance of 9 = /i/7 in Pi [cf. Eqs. (|30| ) -(|32| ) 1 and in f2. Each one of 
the two afforementioned effects is a different reflection of dissipation present in the system. 

A simple case is the self-diffusion problem, i.e., when the tracer and gas particles are 
mechanically equivalent. In this case, = PJ and 7 = 1. This situation has been previously 
studied by Campbell p by means of molecular dynamics simulations. As has been discussed 
in Sec. |HI| , our predictions for the self-diffusion tensor agree qualitatively well with these 
simulations. On the other hand, when the tracer and gas particles are mechanically different, 
to my knowledge no previous studies on the diffusion tensor under shear flow have been 
made. As pointed out in the Introduction, most of the works on granular mixtures M are 
based on the Chapman-Enskog expansion around a local equilibrium state up to the Navier- 
Stokes order, and therefore, they are restricted to the low-dissipation limit in the simple 
shear flow. In addition, they also assume a single temperature to describe the mixture. A 
more careful calculation which takes into account temperature differences has been recently 
made by Garzo and Dufty [T^JT^] . They have obtained explicit expressions for the transport 
coefficients of a granular binary mixture in terms of the restitution coefficients and the 
parameters of the mixture. Since these results have been derived taking the freely cooling 
state as the reference one, the diffusion is characterized by a single scalar coefficient that 
cannot be directly compared to the diffusion tensor obtained here. However, it would be 
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interesting to compare the diffusion results obtained here in the driven sheared case with 
those found in the unforced case [13]. In Fig. [|, we have compared the behavior of the scalar 



\D^ k (which can be understood as a generalized mutual diffusion coefficient in a sheared 
mixture) with the (reduced) diffusion coefficient D* obtained in Ref. [13|] in the tracer limit 



(xi —> 0). We observe that, although the reference states in both descriptions are very 
different, the dependence of both diffusion coefficients on dissipation is quite similar since 
they increase as a increases. This trend is more significant in the unforced case than in the 
sheared case. 

The evaluation of the diffusion tensor for practical purposes requires the truncation of a 
Sonine polynomial expansion. In the case of the (pure) simple shear flow problem, recent 
Monte Carlo simulations |T7| , p3| have shown the accuracy of the leading order truncation. 



We expect that this agreement may be extended to the elements of the diffusion tensor for 
a wide range of values of dissipation. Exceptions to this agreement could be disparate mass 
binary mixtures (e.g., electron-proton systems) for which the first Sonine solution is not 
perhaps a good approximation and higher-order terms should be considered. We hope that 
the results derived in this paper will stimulate the performance of computer simulations 
to check the quality of the approximations given here for the diffusion tensor. Given the 
difficulties associated with molecular dynamics simulations in the low-density regime, one 



could perhaps use the direct simulation Monte Carlo method [26| which is being shown to 
be fruitful in the context of granular systems. 
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APPENDIX A: EVALUATION OF (i, A, B, AND tl 

In this Appendix we evaluate the cooling rates Q, the collisional velocity moments A, 
and B, and the tensorial collision frequency f2 by using the corresponding leading Sonine 
approximations. 

1. Evaluation of Q 

The cooling rate £i of the tracer particles is defined as 

Ci = / dVimxVlMVtlhJi) . (Al) 

A useful identity for an arbitrary function h(Vi) is given by 

Jdv 1 h(v 1 )j 12 [Vi\A°\f2\ =° d i2 1 J dVi J dv 2 fl°\y l )f 2 (y 2 ) 



x J d*Q(& • g)(<x ■ g) [h(V'D - MVO] , (A2) 
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with 



V"i = Vi - ^21(1 + a l2 ){a ■ g)a . 
Using (|A2|) , Eq. ( |70!| ) can be written as 



(A3) 



Ci 



" h r a^ 1 (l + a 12 )fi 21 J dVi J dV 2 f 1 {V 1 )f 2 (V 2 ) J da 6(5- ■ g 12 )(<r ■ g 12 ) 2 



dn{T 

x [/i 21 (l + a 12 )(a ■ g i2 ) - 2(Vi ■ a)] 



To perform the angular integrations, we need the results 

J d*e(a-g 12 )(a-g 12 ) n = p n g? 2 , 



(A4) 



(A5) 



where 



J d&e(& ■ g 12 )(5- ■ g 12 ) n 5- = /3„,+i< 2 



_ (d _ 1)/2 r((n + l)/2) 
Pn r((n + d)/2)" 



n-1. 



>12, 



(A6) 



(A7) 



Thus, the integration over <x in Eq. 



leads to 



Ci 



dn\T] 



r ( 7 1 d 2 - 1 /3 3 (l + a 12 )/i 21 1 dVi | dV 2 h{Y x )f 2 {V 2 ) 



x 



H2igl 2 (l - "12) + 2pi 2 (gi2 • G 



12, 



(A8) 



where Gi 2 = yUi 2 Vi + fi 2 {V 2 is the center of mass velocity. Now, we take the Sonine 



approximations given by (16) for f 2 and (|28|) for f\. Neglecting nonlinear terms in the 
tensors C, = P* — 1, one gets 



where 



7T 



d\l / dV^ 2 /x 21 ^(l-a 12 ) + 2(gt 2 -GI 2 ; 



(A9) 



(A10) 



H ere, V* = Vi/v , g{ 2 = gu/v , G* 12 = G 12 /v , 9 = //i 2 /(// 2 i7) , 7 = 7i/T 2 , and w = 
2T 2 /m 2 . In Eqs. ( |/^9| ) and ( |A10| ), use has been made of the fact that the scalar d cannot 



be coupled to the traceless tensor Q so that the only contributions to Ci come from the 
(pure) Maxwellian terms in (|T^) and (^). The integral 1^ can be evaluated by the change 
of variables 

x = Vt-V 2 *, y = 0V? + V*, (All) 
with the Jacobian (1 + 6)~ d . The integral I({0) can be now easily computed with the result 
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m = r( r(d/2) /2) r(d+3)/2(1 + 9)1/2 [2 ~ 7X21(1 + + M] 



(A12) 



Use of the result dM2|) in Eq. fl£|) yields 

Ci 



2 n 2 af 2 ^0^21(1 + Q12) 



1/2 



[2-/i 21 (l + #)(l + a 12 ) 



The corresponding expression for (2 can be easily obtained from ( |A13| ) and the result is 



(A13) 



n (d-l)/2 



(A14) 



2. Evaluation of A and B 

Since the tensor A can be easily obtained from the expression of B, let us explicitly 
evaluate the latter tensor. It is defined as 



B = y dVimiViViJ^iyil/i,^] 
Using (| A2|) , Eq. ( |A15| ) can be written as 

B = -m 1 a d 2 1 (3 3f i 21 (l + a 12 ) J dV x J d^ 2 h{Y 1 )f 2 {V 2 )g 12 



(A15) 



x 



gi 2 Gi2 + Gi2gi2 + -j^—{2d + 3 - 3tti 2 )gi 2 gi2 - t^7( 1 + a ^)g\ 2 ^ 
d + 3 d + 3 



(A16) 



where use has been made of (|A6|) and 



J d&Q(& ■ g 12 )(S- • g 12 r^ = ^^9l2~ 2 (™gl2gl2 + ^ 2 l) • (A17) 

Substituting the Sonine approximations (|16|) and ( pip for f 2 and /1, respectively, and re- 
taining only linear terms in the tensors Q, one gets 

*„* „-8V-,* 2 -VJ 2 



B = -m 1 a d 2 1 (3 3t x 21 {l + a 12 )n 1 n 2 v 3 o e d/2 7T~ d [ rfV* / dV* 2 g* 12 e~ ev ^ 



x 



1 + eCx-. (v^v* - Jv7 2 u) + c 2 : (v;v; - W 2 2 n 



* /• -1 * /~(* * I /^21 fnj i q q \„* „* A* 21 i \ *2 1 

?12 kj 12 + ^19&12 + 



d + 3 



( 2 rf + 3-3a 12 )g* 2 g* 



12 



d + 3 



1 + a n )g<&$. 



(A18) 



This integral can be computed by the change of variables ( |A11| ). After a lengthy algebra, 
the result is 
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B = — , ,_ (7 12 TOinin 2 )U2i(l + "12)^0 



'1 + 0' 



3/2 



dT(d/2) 
A12 



d 



1 



d + 2 d + 3 
(d+3) 



^21(1 + an) 



11-2- 



2(^ + 2) 



2(d + 2) 



i + 9)\ 12 p; 



(A19) 



where 



A 



12 



1+0 d + 3 



/i 2 i(l + a u ). 



(A20) 



The corresponding expression for A can be easily inferred from Eq. (|A19|) by just making 
the change 1 — > 2 and 9 — *■ 1: 

A = ^ — a{- x m 2 n\vl{\ + a 22 ) {[d + 1 + a 22 (d - 1)]U - (2d + 3 - 3a 22 )P;} . 



2d(d + 2)r(d/2) 



(A21) 



3. Evaluation of ft 



The tensor Q is defined as 



(A22) 



The evaluation of f2 can be made following similar mathematical steps as above. Thus, 
using (|A1|) and the Sonine approximation to / 2 in Eq. ( |A22|) , one gets 



n = ^n- d a d 12 1 f , 21 (l + a 12 )(3 3 n 2 v 3 6 d / 2 J dY{ J dV* 2 g*z\ 2 V\ 



x 



1 + c 2 : (v;v; - ±k* 2 ii 



2 " 2 d v 2 



e -ev«*-v* 2 



2n 2 n~ d fi 21 (l + a 12 )a d 2 1 9 d / 2+1 v p 3 (l + 9)-^ J dx / dyxx(x + y) 



x 



l + (l+0)- 2 C 2 :(y-0x)(y-0x) 



-bx 2 —cy 2 



(A23) 



where 6 = 0(1 + 0) 1 and c = (1 + 0) 1 . This integral is easily performed, with the result 



2tt^ 1 )/ 2 _ /1 + 0' 

n = drm^ n2Vo{1 + ai2) [~o~ ) 

This expression yields Eq. ([55]) given in the text. 



1/2 
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1 







d+21+0 



(A24) 
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FIGURES 



FIG. 1. Plot of the reduced elements of the pressure tensor of the granular gas as functions of 
the restitution coefficient 022 = a in the two-dimensional case. The solid lines correspond to the 
present theory, the dotted ones to that of Ref. [18] and the dashed ones to that of Ref. [19|. 



FIG. 2. Dependence of the temperature ratio 7 = T1/T2 on the restitution coefficient 
022 = ol\2 = a in the two-dimensional case for w = 2 and three different values of the mass 
ratio fj,. The dashed line corresponds to the prediction given by the theory of Jenkins and Mancini 
[21] in the case fj, = 5. 



FIG. 3. Dependence of the diagonal and off-diagonal elements of the reduced self-diffusion 
tensor D* on the restitution coefficient Q22 = 012 = a in the three-dimensional case. 



FIG. 4. Dependence of the scalar (a) ^D kk , the difference (b) D* x — D* y , and the off-diagonal 
elements (c) D* y and (d) D* x of the reduced tracer diffusion tensor D* on the restitution coefficient 
«22 = ot\2 = a in the three-dimensional case for w = 2 and two values of the mass ratio fx: \x = 2 
(solid lines) and fj, = 4 (dashed lines). 



FIG. 5. Plot of the scalar \D kk and the reduced diffusion coefficient D* obtained in Ref. [13] 
as functions of the restitution coefficient 022 = oc\2 = a in the three-dimensional case for w = 2 
and two values of the mass ratio (jl: (jl = 2 (solid lines) and \i = 4 (dashed lines). 
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